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The cosmic microwave background (CMB) is a rich source of cosmological information. Thanks to 
the simplicity and linearity of the theory of cosmological perturbations, observations of the CMB's 
polarization and temperature anisotropy can reveal the parameters which describe the contents, 
structure, and evolution of the cosmos. Temperature anisotropy is necessary but not sufficient to 
fully mine the CMB of its cosmological information as it is plagued with various parameter degenera- 
cies. Fortunately, CMB polarization breaks many of these degeneracies and adds new information 
and increased precision. Of particular interest is the CMB's B-mode polarization which provides 
a handle on several cosmological parameters most notably the tensor-to-scalar ratio, r, and is sen- 
sitive to parameters which govern the growth of large scale structure (LSS) and evolution of the 
gravitational potential. These imprint CMB temperature anisotropy and cause E-to-B-mode po- 
larization conversion via gravitational lensing. However, both primordial gravitational-wave- and 
secondary lensing-induced B-mode signals are very weak and therefore prone to various foregrounds 
and systematics. In this work we use Fisher-matrix-based estimations and apply, for the first time, 
Monte-Carlo Markov Chain (MCMC) simulations to determine the effect of beam systematics on 
the inferred cosmological parameters from five upcoming experiments: PLANCK, POLARBEAR, 
SPIDER, QUIET+CLOVER and CMBPOL. We consider beam systematics which couple the beam 
substructure to the gradient of temperature anisotropy and polarization (differential beamwidth, 
pointing offsets and ellipticity) and beam systematics due to differential beam normalization (differ- 
ential gain) and orientation (beam rotation) of the polarization-sensitive axes (the latter two effects 
are insensitive to the beam substructure). We determine allowable levels of beam systematics for 
given tolerances on the induced parameter errors and check for possible biases in the inferred param- 
eters concomitant with potential increases in the statistical uncertainty. All our results are scaled 
to the 'worst case scenario'. In this case, and for our tolerance levels the beam rotation should not 
exceed the few-degree to sub-degree level, typical ellipticity is required to be 1%, the differential gain 
allowed level is few parts in 10 3 to I0 4 , differential beam width upper limits are of the sub-percent 
level and differential pointing should not exceed the few- to sub-arc sec level. 
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INTRODUCTION 

The standard cosmological model accounts for a mul- 
titude of phenomena occurring over orders of magni- 
tude of length and angular scales throughout the en- 
tire history of cosmological evolution. Remarkably, do- 
ing so only requires about a dozen parameters. Perhaps 
one of the most useful cosmological probes is cosmic 
microwave background (CMB) temperature anisotropy 
whose physics is well understood. Complementary cos- 
mological probes can assist in breaking some of the 
degeneracies inherent in the CMB and further tighten 
the constraints on the inferred cosmological parameters. 
Temperature anisotropy alone cannot capture all the cos- 
mological information in the CMB, and its polarization 
probes new directions in parameter space. B-mode po- 
larization observations are noise-dominated but the ro- 
bust secondary signal associated with gravitational lens- 
ing, which is known up to an uncertainty factor of two 
on all relevant scales, is at the threshold of detection 
by upcoming CMB experiments. The lensing signal may 



have been detected already through its signature on the 
CMB anisotropy as reported recently by ACBAR (Re- 
ichardt et al. [1]). Lensing by the large scale structure 
(LS) also converts primordial E-mode to secondary B- 
mode. When high fidelity B-mode data are available a 
wealth of information from the inflationary era (Zaldar- 
riaga & Seljak [2], Kamionkowski, Kosowsky & Stebbins 
[3]), and cosmological parameters that control the evo- 
lution of small scale density perturbations (such as the 
running of the spectral index of primordial density per- 
turbations, neutrino mass and dark energy equation of 
state), will be extracted from the CMB. At best, B-mode 
polarization from lensing is a factor of three times smaller 
than the primordial E-mode polarization, so it is prone 
to contamination by both astrophysical foregrounds and 
instrumental systematics. It is mandatory to account for, 
and remove when possible, all sources of spurious B-mode 
in analyzing upcoming CMB data, especially those gen- 
erated by temperature leakage due to beam mismatch, 
since temperature anisotropy is several orders of magni- 
tude larger than the expected B-mode level produced by 
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lensing. 

Beam systematic have been discussed extensively (Hu, 
Hedman & Zaldarriaga [4], Rosset et al. [5], O'Dea, 
Challinor & Johnson [6], Shimon et al. [7]). All the 
effects are associated with beam imperfections or beam 
mismatch in dual beam experiments, i.e. where the po- 
larization is obtained by differencing two signals which 
are measured simultaneously by two beams with two or- 
thogonal polarization axes. Fortunately, several of these 
effects (e.g. differential gain, differential beam width 
and the first order pointing error - 'dipole'; Hu, Hed- 
man & Zaldarriaga [4], O'Dea, Challinor & Johnson [6], 
Shimon et al. [7]) are reducible with an ideal scanning 
strategy and otherwise can be cleaned from the data set 
by virtue of their non-quadrupole nature which distin- 
guishes them from genuine CMB polarization signals. 
Other spurious polarization signals, such as those due 
to differential ellipticity of the beam, second order point- 
ing errors and differential rotation, persist even in the 
case of ideal scanning strategy and perfectly mimic CMB 
polarization. These represent the minimal spurious B- 
mode signal, residuals which will plague every polariza- 
tion experiment. We refer to them in the following as 
'irreducible beam systematics'. We assume throughout 
that beam parameters are spatially constant. Two re- 
cent works (Kamionkowki [8] and Su, Yadav & Zaldar- 
riaga [9]) considered the effect of spatially-dependent sys- 
tematic beam-rotation and differential gain, respectively. 
This scale-dependence and the associated new angular 
scale induce non-trivial higher order correlation functions 
through non-gaussianities which can be both used to op- 
timally remove the space-dependent component of beam 
rotation [8] and mimic the CMB lensing signal, thereby 
biasing the quadratic estimator of the lensing potential 
[9]. 

To calculate the effect of beam systematics we invoke 
the Fisher information-matrix formalism as well as Monte 
Carlo simulations of parameter extraction, the latter for 
the first time. Our objective is to determine the suscep- 
tibility of the above mentioned, and other, cosmological 
parameters to beam systematics. For the Fisher-matrix- 
based method and the Monte Carlo simulations we calcu- 
late the underlying power spectrum using CAMB (Lewis, 
Challinor & Lascnby [10]). The Monte-Carlo simulations 
are carried out with COSMOMC (Lewis & Bridle [11]). 
We represent the extra noise due to beam systematics by 
analytic approximations (Shimon et al. [7]) and include 
lensing extraction in the parameter inference process, fol- 
lowing Kaplinghat, Knox & Song [12] and Lesgourgues et 
al. [13] (see also Perotto et al. [14] for the Monte Carlo 
simulations) for neutrino mass (and other cosmological 
parameters) reconstruction from CMB data. 

This paper illustrates the effect of beam systematics 
and its propagation to parameter estimation and error 
forecasts for upcoming experiments. Our main concern 
is the effect on the following cosmological parameters: 



the tensor-to-scalar ratio r, the total neutrino mass M v 
(assuming three degenerate species), tilt of the scalar in- 
dex a, dark energy equation of state w, and the spatial 
curvature, Qk- The lensing- induced B-mode signal is sen- 
sitive to all parameters (except the tensor-to-scalar ratio) 
and peaks at few arcminute scales, while the tensor-to- 
scalar ratio depends on the energy scale of inflation and 
the primordial signal peaks at the characteristic horizon 
size at last scattering, w 2°. We note that while the LSS- 
induced and primordial tensor power B-mode spectra are 
sub-^iK the shape of the primordial B-mode spectrum is 
known (only its amplitude is unknown, Keating [15]) and 
the secondary LSS-induced B-mode is guaranteed to ex- 
ist by virtue of the known existence of LSS and E-mode 
polarization. 

The paper is organized as follows. We describe the 
formalism of beam systematics for general non-gaussian 
beams and provide a cursory description of a critical tool 
to mitigate polarization systematics -a half wave plate 
(HWP), in section 2. The effect of lensing on parame- 
ter extraction within the standard quadratic-estimators 
formalism is discussed in section 3. The essentials of the 
Fisher matrix formalism are given in section 4 as well 
as some details on the Monte Carlo simulations invoked 
here. Our results are described in section 5 and we con- 
clude with a discussion of our main findings in section 
6. 



BEAM SYSTEMATICS 

Beam systematics due to optical imperfections depend 
on both the underlying sky, the properties of the po- 
larimeter and on the scanning strategy. Temperature 
anisotropy leaks to polarization when the output of two 
slightly different beams with orthogonal polarization- 
sensitive directions is being differenced. A trivial exam- 
ple is the effect of differential gain. If the two beams have 
the same shape, width, etc. except for different overall 
response, i.e. normalization, the difference of the mea- 
sured intensity will result in a non- vanishing polarization 
signal. Similarly, if two circular beams slightly differ by 
their width this will again induce a non- vanishing polar- 
ization upon taking the difference (see Fig. 2 at Shimon 
et al. [7]). The spurious polarization will be propor- 
tional to temperature fluctuations on scales comparable 
to the difference in beamwidths, which, due to the cir- 
cular symmetry of the problem, will be proportional to 
second order gradients of the temperature anisotropy. To 
eliminate these effects this beam imperfection has to cou- 
ple to non-ideal scanning strategy as described in Shimon 
et al. [7] and below. A closely related effect, which does 
not couple to scanning strategy, is the effect of differ- 
ential beam ellipticity. Here, the spurious polarization 
scales as the second order gradient of the temperature 
anisotropy to leading order. Another effect, widely de- 
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scribed in the literature, which again couples beam asym- 
metry and temperature anisotropy to scanning strategy 
is the effect due to differential pointing. The idea is sim- 
ple; if two beams point at two slightly different directions 
they will statistically measure two different intensities 
proportional to fluctuations of the background radiation 
on these particular scales. The difference, which may be 
naively regarded as polarization, is non- vanishing in this 
case provided the scanning strategy is non-ideal and con- 
tains either a dipole and/or an octupole (Shimon et al. 
[7]). Finally, the effect of beam rotation we consider in 
this work is due to uncertainty in the overall beam orien- 
tation. This mixes the Q and U Stokes parameters and as 
a result also leaks E to B and vice versa. A constructive 
order-of-magnitude example is the effect of differential 
pointing. This effect depends on the temperature gradi- 
ent to first order. The rms CMB temperature gradients 
at the 1°, 30', 10', 5' and 1' scales are w 1.4, 1.5, 3.5, 2.5 
and 0.2 /iX/arcmin, respectively. Therefore, any temper- 
ature difference measured with a dual-beam experiment 
(with typical beamwidth few arcminutes) with awl' 
pointing error will result in a w 1/uK systematic polar- 
ization which has the potential to overwhelm the B-mode 
signals. 

Similarly, the systematic induced by differential ellip- 
ticity results from the variation of the underlying tem- 
perature anisotropy along the two polarization-sensitive 
directions which, in general, differ in scale depending on 
the mean beamwidth, degree of ellipticity and the tilt 
of the polarization-sensitive direction with respect to the 
ellipse's principal axes. For example, the temperature 
difference measured along the major and minor axes of 
a 1° beam with a 2% ellipticity scales as the second gra- 
dient of the underlying temperature which on this scale 
is fts 0.2^if/arcmin 2 and the associated induced polar- 
ization is therefore expected to be on the w pK level. 
The spurious signals due to pointing error, differential 
beamwidth and beam ellipticity all peak at angular scales 
comparable to the beam size. If the beam size is « 1° 
the beam systematics mainly affect the deduced tensor- 
to-scalar ratio, r. If the polarimeter's beamwidth is a 
few arcminutes the associated systematics will impact the 
measured neutrino mass m„, spatial curvature fife, run- 
ning of the scalar spectral index a and the dark energy 
equation of state w (which strongly affects the lensing- 
induced B-mode signal). It can certainly be the case that 
other cosmological parameters will be affected as well. 

Two other spurious polarization signals we explore are 
due to differential gain and differential rotation; these ef- 
fects are associated with different beam 'normalizations' 
and orientation, respectively, and are independent of the 
coupling between beam substructure and the underly- 
ing temperature perturbations. In particular, they have 
the same scale dependence as the primordial temperature 
anisotropy and polarization power spectra (as long as the 
differential gain and beam rotation are spatially inde- 



pendent; this of course changes if they depend on space 
[8], [9]), respectively, and their peak impact will be on 
scales associated with the CMB's temperature anisotropy 
(~ 1°) and polarization (« 10'). 

Mathematical Formalism 

We work entirely in Fourier space and in this section 
we generalize our results (Shimon et al. [7]) to the case of 
the most general beam shapes. Although the tolerance 
levels on the beam parameters we derive in sections 4 
and 5 are based on the assumption of elliptical beams, 
they can be easily generalized to arbitrary beam shape, 
given the beam profile, as we describe below. This can be 
used to adapt our results to actual measured beam maps 
incorporating other classes of beam non-ideality such as 
sidelobes. 

We expand the temperature anisotropy and Q and U 
Stokes parameters in 2-D plane waves since for sub-beam 
scales this is a good approximation. While the (spin- 
0) temperature anisotropy is expanded in scalar plane 
waves e ll ' r , the (spin ±2) polarization tensor Q + ill is 
expanded in tensor plane waves e A- r e ± 2l (4>\-4>r) wnere <^ r 
is the angle defining the direction of the radius-vector r 
in real space as conventional (in an arbitrarily coordinate 
system on the sky (j> r is the azimuthal angle along the 
line of sight) and cf)\ defines the direction of the wave- 
vector 1 in 1-space in a coordinate-system fixed to the 
beam as defined below, in Eq. (2). Since in real space the 
temperature and polarization fields are convolved with 
the polarimeter's beams, these expressions are simply the 
product of their Fourier transforms in Fourier space. For 
a general beam B{v) the measured 2-D polarized beam 
map may exhibit a pointing error p. In this case, the 
Fourier transform of the beam function acquires a phase 

B(l)-B(l)exp(il-p). (1) 

It is useful to switch to polar coordinates at this point, 
where we define 

l x = I cos((/>i + ip — a) 
l y = I sin(0i + ip — a) 
p x = peon 9 

p y = p sin 6 (2) 

and a = f3 + 8 + ip is the angle of the polarization axis in 
a coordinate system fixed to the sky (Fig. 1 of Shimon et 
al. [7] ) . The Fourier representation of an arbitrary beam 
then becomes 

B(l) = J B(r)e a r d 2 r 

_ J _g( r )g iir cos(0 ; +-t/J-a-</v)+i(pcos(0 ; +0+i/;-a-<6r) 

= Y. B ^ i y (m+n){ ^ a) ( 3 ) 

m.n 
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where 

x J rdrJ rn {lr) j d<f) r B{v)e^ m+n ^ , (4) 

and in the last step we employed the expansion of 2-D 
plane waves in terms of cylindrical Bessel functions 

n—oo 

e il P co S (4>i-<t> P ) = ^ i n J n (lp)e m{4 ''- 4 ' p) . (5) 

n— — oo 

As in Shimon ct al. [7], the optimal map constructed 
from the CMB data depends on the measurements as 

m(p) = Aj(p,t)A 1 (p,t) 
\t,jep 



Aj(p,t)d J (p,t) 



(6) 



where the sums run over all measurements of the pixel 
p. The pointing vector A is given by 



1 

'2 l 



1 

'2' 



-2ia(p,t] 



(7) 



a is a function of both the pixel p and t, and A T is A 
transposed. Once the leading beam coefficients B mn (l) 
have been calculated, the induced power spectra of the 
systematics can be calculated according to Eqs. (24), 
(33), (A.l) and (A.2) of Shimon et al. [7]. 

Several of the beam systematics can be mitigated by 
employing a rotating half wave plate (HWP) polariza- 
tion modulator (e.g. Hanany et al. [16], Johnson et al. 
[17], MacTavish et al. [18]). These can operate in con- 
tinuous or stepped rotation. When HWP modulators are 
included we replace the above scanning angle a(p, t) with 
a(p, t) + 2ipt where ip is the angular velocity of the HWP 
(O'Dea, Challinor & Johnson [6]). Our deduced tolerance 
levels given below are presented in a fashion independent 
of the details of the scanning strategy; all the informa- 
tion about the scanning strategy is encapsulated in the 
functions /i, / 2 and f?,: 



fi 
h 

h 



^IM-i,o)| 2 

l|M-i,-i)l 2 + ^IM-M)l 2 
\ (/(o,i)^(i,-i)) 



(8) 



where 



/(m,n) 
h±(m, n) 
D 



| e -t(2ra+n)a(p,t)| 

1 [/(™,n)-/(m±2,n)(e ±4M ( p ^)] 



D 

1 - (e 



and the angular brackets in (e mQ ( p '*)) represent average 
over measurements of a single pixel p, averaged over time. 
In these averages a(p, t) — ► a(p,i) + 2<pt, and therefore 
even if the scanning strategy does not uniformly cover all 
polarization angles a of a given spatial pixel, the HWP 
mitigates the spurious polarization caused by beam sys- 
tematic effects if integrated over long time intervals. If 
the hexadecapole of the scanning strategy is negligible, 
the scanning strategy function f\ depends only on the 
quadrupole moment of the scanning strategy while / 2 en- 
capsulates information on both the dipole and octupolc 
moments of the scanning strategy. 

At this point, it is instructive to show how, in the case 
of an ideal scanning strategy the first order pointing ef- 
fect vanishes. As can be seen from Table II this effect 
involves a convolution of the beam function and under- 
lying temperature anisotropy power spectrum with / 2 in 
multipole space. From the above definitions 



M-i, -i) 
M-M) 
Mi,-i) 



1/D[(e 3la ) 
1/D [(e ia ) (e 

/*;(-!, i) 



e- ia )(e 4ia )] 
3ia )(e 4ta )] 



(10) 



-i4a(p,t)^ e i4a;(p,t) 



(9) 



and therefore if with each scanning angle a there is asso- 
ciated an angle a + 180° the / 2 (Eqs. 8 and 10) vanishes 
in real space and so does its Fourier transform. Note that 
even if the scanning strategy is non-ideal / 2 will vanish 
provided that for each angle a the angle a + 180° is sam- 
pled the same number of times per pixel. This suggests 
that the dipole systematic can be completely removed by 
removing all data points that contribute to /i+(— 1, — 1) 
and h + (— 1,1), i.e. those measurements at a for which 
a + 180° is not sampled. Similar considerations apply to 
fi which controls the level of the differential beamwidth- 
and differential gain-induced systematics (see Table II). 



Simplifying Scan Strategy Effects 

When the polarization angle at each pixel on the sky is 
uniformly sampled the average (e ma ) vanishes for every 
n^0. In this case the scanning strategy is referred to as 
an ideal scanning strategy. For uniform, but non-ideal, 
scanning strategies, the scanning functions fi, f 2 and 
/ 3 mentioned above (which are combinations of (e mQ )) 
are non-vanishing even when n > but uniform in real 
space. As a result their Fourier transforms are unnormal- 
ized delta-functions (the actual amplitudes are directly 
related to the average values (e ma )), and the convolu- 
tions in Fourier space shown in Tables III-IV of Shimon 
et al. [7] become trivial. To determine the tolerance level 
for beam parameters we assume such uniform scanning 
strategies. 

A uniform scanning strategy is a particularly useful ex- 
ample. A nearly-uniform scanning strategy can be rea- 
sonably approximated by a sum over a few lowest multi- 
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poles, such as 



fi{i) = Y,fs'w-i')/{i-i'). 



(ii) 



J'=0 



Here \l ma x\ is assumed sufficiently small, and i = 1 or 2, 
where /i, /2 and / 3 are defined in Eq. 8). In this case the 
£-mode mixing due to the convolution of the underlying 
power spectra and the scanning functions as in Table II 
(* stands for convolution in multipole space), fi can be 
written as 



A 

2tt 



■C 



(12) 



We assume here that the nonvanishing multipoles of /j 
are concentrated near 0, i.e. that the scanning strategy is 
non-ideal, yet approximately uniform. We have employed 
this simplifying assumption throughout. 



THE EFFECT OF SYSTEM ATICS ON LENSING 
RECONSTRUCTION 

Gravitational lensing of the CMB is both a nuisance 
and a valuable cosmological tool (e.g. Zaldarriaga & Sel- 
jak [19]). It certainly has the potential to complicate 
CMB data analysis due to the non-gaussianity it induces. 
However, it is also a unique probe of the growth of struc- 
ture in the linear, and mildly non-linear, regimes (rcdshift 
of a few). Kaplinghat, Knox & Song [12], Lesgourgues et 
al. [13], as well as others, have shown that with a nearly 
ideal CMB experiment (in the sense that instrumental 
noise as well as astrophysical foregrounds are negligibly 
small), neutrino mass limits can be improved by a factor 
of approximately four by including lensing extraction in 
the data analysis using CMB data alone. This lensing 
extraction process is not perfect; a fundamental residual 
noise will afflict any experiment, even ideal ones. This 
noise will, in principle, propagate to the inferred cosmo- 
logical parameters if the latter significantly depend on 
lensing extraction, e.g. neutrino mass, a and w. It is 
important to illustrate first the effect of beam system- 
atics on lensing reconstruction. By optimally filtering 
the temperature and polarization Hu & Okamoto [20] re- 
constructed the lensing potential from quadratic estima- 
tors. It was shown that for experiments with ten times 
higher sensitivity than Planck, the EB estimator yields 
the tightest limits on the lensing potential. This result 
assumes no beam systematics which might significantly 
contaminate the observed B-modc. 

We illustrate the effect of differential beam rotation, 
ellipticity and differential pointing (see Shimon et al. 



[7]) on the noise of lensing reconstruction with POLAR- 
BEAR (1200 detectors), CMBPOL-A (one of two toy 
experiments we consider for CMBPOL; 0. 22yU.fr sensi- 
tivity and 5' beam) and a toy- model considered ear- 
lier by O'Dea, Challinor & Johnson [6] which we refer 
to as QUIET+CLOVER in Figures 1, 2 and 3, respec- 
tively. These are perhaps the most pernicious system- 
atics. Beam rotation induces cross-polarization which 
leaks the much larger E-mode to B-mode polarization 
and differential ellipticity leaks T to B. Both leak to B- 
mode in a way indistinguishable under rotation from a 
true B-mode signal. The rotation and ellipticity param- 
eters (e and e, respectively) we considered range from 
0.01 to 0.20 (e is dimensionless and e is given in radi- 
ans). The differential pointing p, was set to 1% and 10% 
of the beamwidth while the dipole and octupole compo- 
nents of the scanning strategy were set to the 'worst case 
scenario' = 2ir, i.e. the unlikely situation where all 
'hits' at a given pixel take place at the same polarization 
angle a (again, for ideal scanning strategy / 2 = and 
the dipole effect due to differential pointing vanishes). 
Note for POLARBEAR (Figure 1) with e,e = 0.2 the 
lensing potential can be reconstructed up to I « 200, 
while with no beam rotation it can be reconstructed up to 
I « 250. However, with CMBPOL-A (Figure 2) lensing 
reconstruction degrades significantly in the presence of 
beam rotation (from good reconstruction up to 1«600 in 
the systematics-free case down to 1« 250 when e, e = 0.2 
and p = 0.5' (in case f 2 = 2tt)). The reason for the qual- 
itative difference is that for experiments with sensitivi- 
ties comparable to PLANCK or POLARBEAR, the best 
estimator of the lensing potential comes from the TT, 
TE and EE correlations (depending on scale I) and the 
cross-correlations involving B-mode are only secondary in 
probative power (see top left panel of Fig.l). Therefore, 
lensing reconstruction for these experiments is hardly af- 
fected by beam systematics (we ignored the negligible 
beam systematics' effect on temperature anisotropy and 
considered only those of E and B). In contrast, as can be 
seen from Figure 2, CMBPOL-A's lensing reconstruction 
is significantly degraded since its lensing reconstruction 
is dominated by the contribution of the EB estimator for 
all relevant multipoles (see top left panel of Fig. 2). The 
modified noise in reconstructing the lensing potential, 
Nf- d , is consistently substituted into our Fisher matrix 
and Monte Carlo simulations (below, we summarize the 
relevant expressions of the quadratic estimators method) . 

Following Hu & Okamoto [20] the MV noise on the 
lensing deflection angle reconstructed power spectrum 
Cf d is 



AT 



dd 
MV 



)a(3 



a/3 



(13) 
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where 



N af3 (L) = L- 2 A a {L)A p {L) J A^F a (li,l 2 ) (F^(li,l 2 ) 



h h 



(27T) 



_X Q ,x, _ X ' )X , 



A a (L) = i 2 



(2^ 



n -1 



^(11,11)^(11,11) 



(14) 



and a stands for one of the pairings TT,TE,EE,TB and 
EB (BB does not participate in these combinations). The 
coupling takes place between different modes l\ and 1%. 
When a = TT or EE 



and when a = TB or EB 



h 1-2 



rjixxQix 1 x' 



(15) 



(16) 



where C[ are the observed power spectra, i.e. including 
lensing, main-beam dilution on small scales, and in prin- 
ciple - beam systematics (see Tables I and II). The latter 
mainly affect the B-mode polarization, and as a result, 
the EB estimator. A list of (li , I2) can be found in 
Hu & Okamoto [20]. We have used the publically avail- 
able code [21] employed in Lesgourgues et al. [13] and in 
Perotto et al. [14]. The code is based on the formalism 
developed in Okamoto & Hu [22], an extension of Hu & 
Okamoto [20] to the full-sky, to calculate the noise level 
in lensing reconstruction. 



ERROR FORECAST 

Accounting for beam systematics in both Stokes pa- 
rameters and lensing power spectra is straightforward. In 
addition, the instrumental noise associated with the main 
beam is accounted for, as is conventional, by adding an 
exponential noise term. Assuming gaussian white noise 



N, 



1 



(17) 



where a runs over the experiment's frequency bands. The 
noise in channel a is (assuming a gaussian beam) 



iVr = (^A a ) 2 e /(i+1)e ' /81n(2) , 



(18) 



where A a is the noise per pixel in /uK-arcmin, 9 a is the 
beam width (see Table III), and we assume noise from 
different channels is uncorrelated. The power spectrum 
then becomes 



C, 



x 



C, 



x 



N 



(19) 



where X is either the auto-correlations TT, EE and BB 
or the cross- correlations TE, TB and EB (the latter two 



power spectra vanish in the standard model but not in the 
presence of beam systematics and exotic parity- violating 
physics e.g. Carroll [23], Liu, Lee & Ng [24], Xia et al. 
[25], Komatsu et al. [26] or primordial magnetic fields 
e.g. Kosowsky & Loeb [27]). For the cross-correlations, 
the N x vanish as there is no correlation between the 
instrumental noise of the temperature and polarization 
(in the absence of beam systematics). 



Fisher-Matrix-Based Calculation 

The effect of instrumental noise is simply to increase 
the error bars, which is evident from the Fisher matrix 
formalism below. The l-o~ error o~(\i) on the cosmolog- 
ical parameter A^ can be read-off from the appropriate 
diagonal element of the inverse Fisher matrix 



where the Fisher matrix elements are defined as 

d 2 L 



F ■■ — — 



(20) 



(21) 



L is the likelihood function, and Eq. (|2ip is evaluated 
at the best-fit point in parameter space. Explicitly, the 
Fisher matrix elements for the CMB read 



F H = |^(2^ + l)/ sky Trace 



c-xacc-xflc 



(22) 



The pre-factor ^(21 + l)/ s ky comes from the sample- 
variance of the multipole I with an experiment covering 
a fraction / s k y of the sky. The matrix C is 



/ C' TT 





'I'd 
I 



qiEE 









(jiBB 






Qidd 



\ 



(23) 



where the diagonal primed elements C[ xx = 
N xx and X € {T, E, B, d}. In general, N t BE = N t BB = 
2Nf T ■ Note that, except for Nf d , which is not an in- 
strumental noise and emerges only because of the limited 
reconstruction of the lensing potential by the quadratic 
estimators of Hu & Okamoto [20] , the instrumental noise 
will increase C, but not its derivatives with respect to 
the cosmological parameters. This will increase the er- 
ror on the parameter estimation as seen from Eqs. (|20[) . 
(f2"2")l . and (f2"3"]) . It is merely because the instrumental 
noise dilutes the information below the characteristic 
beamwidth scale, and the error increases correspondingly. 
However, this is not necessarily the case with beam sys- 
tematics since they couple to the underlying cosmological 
model, and therefore do depend on cosmological param- 
eters. This noise due to systematics, iV ; sys , contributes 
to both C and and its effect on the confidence level 
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FIG. 1: For all panels the solid black curve is the deflection angle power spectrum C; caused by gravitational lensing by the 
LSS. Top left: The noise (with no systematics) in lensing-reconstruction from the quadratic optimal filters for POLARBEAR TT 
(solid dark blue), EE (dot light blue), TE (dashed green), TB (dot-dash yellow), EB (double-dot-dash orange) and MV (dashed 
dark red) . For POLARBEAR sensitivity and angular resolution, the lowest-noise estimator is one of the EE and TE estimators 
depending on the angular size. Top right: Noise in lensing reconstruction for POLARBEAR with the MV estimator including 
the effects of the most pernicious irreducible cross-polarization systematic: differential rotation. Differential rotation values 
are (bottom to top) e =0.01, 0.02, 0.05, 0.10 and 0.20 radian, respectively. High signal-to-noise deflection angle reconstruction 
can be obtained over nearly a decade of angular scale. The lensing reconstruction is not significantly affected by systematics 
because of the significant contribution of the temperature to the MV estimator. Bottom left: Noise in lensing reconstruction 
for POLARBEAR with the MV estimator including the most pernicious irreducible instrumental polarization systematics: 
differential ellipticity. Differential ellipticity values are (bottom to top) e =0.01, 0.02, 0.05, 0.10 and 0.20, respectively (we 
assume ip = 45°). The same explanation for insensitivity to differential rotation applies here for differential ellipticity e; the 
best estimator for this experiment is derived from temperature correlations which are hardly affected by beam systematics (and 
completely ignored in this analysis). Bottom right: Noise in lensing reconstruction for POLARBEAR with the MV estimator 
for the most pernicious reducible instrumental polarization systematics: differential pointing with 1% and 10% pointing errors 
(i.e. p — 0.01(7 and p = O.lcr, respectively) under the 'worst case' assumption that the scanning-strategy-related function 
/ 2 = 2tt. 



of parameter estimation can be, in principle, either a 
degradation or an improvement. This argument ignores 
potential systematic errors, i.e., bias (systematic shift of 
the average of the statistical distribution which charac- 
terizes certain cosmological parameters) of the recovered 
average values of the cosmological parameters. Indeed, 
as we show below, the main effect of beam systemat- 
ics is to bias the inferred cosmological parameters, es- 



pecially for large beam mismatch parameters, as naively 
expected (see Fig. 4 for a comparison between the bias 
and the uncertainty induced on the tensor-to-scalar ratio 
and neutrino masses by beam ellipticity). It is important 
to note that, although our focus is on beam systemat- 
ics and their effect on parameter estimation, we do not 
include the systematics-induced Cj B and Cf B in the 
analysis (Eqs. 22, 23) because our main concern is how 
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FIG. 2: Lensing reconstruction with CMBPOL-A: As in Figure 1. For CMBPOL-A sensitivity and angular resolution, the 
lowest-noise estimator comes from correlations of the EB estimator. Therefore, lensing reconstruction is only mildly affected 
by beam systematics. 



standard data analysis pipelines will be affected by beam 
systematics. We defer the treatment of the more gen- 
eral case, which includes the parity-violating terms and 
their effect on constraints of beyond-the-standard-model 
parity-violating interactions in the primordial universe, 
to a future work (Shimon, Miller & Keating [28]). 

Given the beam systematics the bias of a parameter 
can be calculated with the Fisher Matrix. This has been 
done by O'Dea, Challinor & Johnson [6]. The bias in a 
parameter A^, if not too large, is given by 

AX t = (\f s ) - (\m = ]T (F- 1 ) .. B 3 (24) 

3 

where the bias vector B can be written as 

B^tCD'S" 1 ^- (25) 
i j 

and Sy = cov(C;,C/) and C; is a vector containing all 
six power spectra. 



Monte Carlo Simulations 

The Fisher matrix approach is known to provide reli- 
able approximation to the uncertainty in case of gaussian 
distributions and only a lower bound for more general dis- 
tributions by virtue of the Cramer-Rao theorem. It can 
yield poor estimates, however, in cases of large biases and 
parameter degeneracies. To check for such effects in our 
simulations we repeated the analysis with MCMC simu- 
lations which make use of the full likelihood function and 
not only its peak value and can therefore provide reliable 
estimates of parameter errors even in the presence of large 
biases. Our simulations illustrate that even when we con- 
sider the Fisher-matrix-based results as a guide for choos- 
ing the beam parameters for the MCMC simulations, the 
resulting bias that the Monte-Carlo simulations predict 
can be larger than those found with the Fisher-matrix- 
based calculation (actually Eq. 24 assumes that the bias 
is small compared to the characteristic width of the like- 
lihood function of the parameter in question; when this is 
not the case this approximation is invalid) in some cases. 
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FIG. 3: Lensing reconstruction with QUIET+CLOVER: As in Fig.l. For QUIET+CLOVER sensitivity and angular resolution, 
the lowest-noise estimator comes from correlations of the EB estimator. Therefore, lensing reconstruction is significantly affected 
by beam systematics. 



This important point is further elucidated in the next sec- 
tion. For our Monte Carlo simulations we use a modified 
version of CosmoMC [29] which includes measurements 
of the lensing potential and its cross-correlation with the 
temperature anisotropy when calculating the likelihood 
in order to run these simulations. An eleven parameter 
model is used (flbh 2 , Qd m h 2 , 6, r, Q„h 2 , w, n s , n tl a, 
log(10 10 A s ), and r). We ran simulations for each of the 
five systematic effects with noise corresponding to PO- 
LARBEAR, CMBPOL-B and QUIET+CLOVER exper- 
iments. While running Monte Carlo simulations is much 
more time consuming compared to the Fisher matrix ap- 
proach, the error forecasts for future experiments will 
ultimately have to account for a potentially significant 
biases of the inferred cosmological parameters. 

RESULTS 

For both the Fisher and MCMC methods we consider 
the effect of both irreducible and reducible systematics. 



By 'reducible' we refer to systematics which depend on 
the coupling of an imperfect scanning strategy to the 
beam mismatch parameters. These can, in principle, be 
removed or reduced during data analysis. This includes 
the differential gain, differential beamwidth and first or- 
der pointing error beam systematics. By 'irreducible' we 
refer to those systematics that depend only on the beam 
mismatch parameters (to leading order). For instance, 
the differential ellipticity and second order pointing er- 
ror persist even if the scanning strategy is ideal. For 
reducible systematics the scanning strategy is a free pa- 
rameter in our analysis (under the assumption it is non- 
ideal, yet uniform, over the map) and we set limits on the 
product of the scanning strategy (encapsulated by the /i 
and parameters) and the differential gain, beamwidth 
and pointing, as will be described below. 

To calculate the power spectra we assume the con- 
cordance cosmological model throughout; the baryon, 
cold dark mater, and neutrino physical energy densities 
in critical density units \7l h 2 = 0.021, fl c h 2 = 0.111, 
Vl v h 2 = 0.006. The latter is equivalent to a total neu- 
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FIG. 4: Uncertainty in the tensor-to-scalar ratio r (top) and 
total neutrino mass M v (bottom) due to beam ellipticity of 
POLARBEAR. The black solid curve is the statistical error 
(uncertainty) and the red dashed curve is the bias. As we 
vary e, the uncertainty increases by only a few percent (i.e. 
the width of the corresponding 1-D likelihood function does 
not significantly change). The bias, however, sharply rises 
with increasing ellipticity, i.e. the expectation value of r and 
M v significantly changes. In general, we find that beam sys- 
tematics mainly bias the inferred parameters since, for large 
enough beam systematics parameters, the spurious polariza- 
tion signal overwhelms the cosmological signal. 



trino mass M v — m v,i =0.56eV, slightly lower than 

the current limit set by a joint analysis of the WMAP 
data and a variety of other cosmological probes (0.66eV, 
e.g. Spergel et al. [30]). We assume degenerate neutrino 
masses, i.e. all neutrinos have the same mass, 0.19 eV, 
for the purpose of illustration, and we do not attempt 
to address here the question of what tolerance levels are 
required to determine the neutrino hierarchy. As was 
shown by Lesgourgues et al. [13], the prospects for deter- 
mining the neutrino hierarchy from the CMB alone, even 
in the absence of systematics, are not very promising. 
This conclusion may change when other probes, e.g. Ly- 
a forest, are added to the analysis. Dark energy makes 
up the rest of the energy required for closure density. The 
Hubble constant, dark energy equation of state and he- 



lium fraction are, respectively, Hq = 70 km sec _1 Mpc _1 , 
W = — 1 and Yffe = 0.24. h is the Hubble constant in 100 
km/sec/Mpc units. The optical depth to reionization 
and its redshift are r re = 0.073 and z re = 12. The nor- 
malization of the primordial power spectrum was set to 
A s = 2.4 x 10~ 9 and its power law index is n s — 0.947 
(Komatsu et al. [26]). 

Since the effect of beam systematics is the focus of 
this paper, and because these systematics are generally 
manifested on scales smaller than the beamwidth (except 
for the effects of differential gain and rotation) we con- 
centrate on some specific cosmological parameters which 
will be targeted by upcoming CMB experiments. These 
parameters have been chosen for our analysis because we 
expect them to benefit from lensing extraction or simply 
because they are somehow associated with small angu- 
lar scales and therefore are prone to systematics on sub- 
beam scales. We limit our analysis to the tensor-to-scalar 
ratio r, dark energy equation of state w, spatial curvature 
flk, running of the scalar index a and total neutrino mass 
M v . While r is mainly constrained by the primordial B- 
mode signal that peaks on degree scales (and is therefore 
not expected to be overwhelmed by the beam systemat- 
ics which peak at sub-beam scales), it is still susceptible 
to the tail of these systematics, extending all the way 
to degree scales, because of its expected small amplitude 
(less than 0.1/iK). The tensor-to-scalar ratio is also af- 
fected by differential gain and rotation which are simply 
rescalings of temperature anisotropy and E-mode polar- 
ization power spectra, respectively, and therefore do not 
necessarily peak at scales beyond the primordial signal. 

The other four parameters either determine the pri- 
mordial power spectrum P(k) on small angular scales 
(e.g. a) or affect the lensed signal (both temperature 
and polarization) at late times (e.g. M v , fife, w). Ideally, 
the lensing signal, which peaks at I « 1000, provides a 
useful handle on the neutrino mass as well as other cos- 
mological parameters which govern the evolution of the 
large scale structure and gravitational potentials. How- 
ever, the inherent noise in the lensing reconstruction pro- 
cess (Hu & Okamoto [20]) which depends, among others, 
on the instrument specifications (instrumental noise and 
beamwidth), now depends on beam systematics as well. 
The systematics, however, depend on the cosmological 
parameters through temperature leakage to polarization, 
and as a result there is a complicated interplay between 
these signals and the information they provide on cosmo- 
logical parameters. As our numerical calculations show, 
the effect on the inferred cosmological parameters stems 
from both the direct effect of the systematics on the pa- 
rameters (the top-left 3x3 block of the covariance matrix, 
Eq. 23) and the indirect effect on the noise in the lensing 
reconstruction, 7V ; dd , in cases where the MV estimator is 
dominated by the EB correlations (see section 3). 
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Fisher Matrix Results 

The Fisher information-matrix gives a first order ap- 
proximation to the lower bounds on errors inferred for 
these parameters. However, by construction, it uses only 
the information from the peak of the likelihood function. 
Markov Chain Monte Carlo simulations are known to 
be superior to Fisher-matrix-based analysis in cases of 
strong parameter degeneracies and bias but Fisher ma- 
trix results are useful for first order approximation and 
provide starting values for MCMC simulations. 

We follow O'Dea, Challinor & Johnson [6] in quantify- 
ing the required tolerance on the differential gain, differ- 
ential beamwidth, pointing, ellipticity and rotation. To 
estimate the effect of systematics and to set the system- 
atics to a given tolerance limit one has to compare the 
systematics-free la error in the i-th parameter (Eq.20) 
to the error obtained in the presence of systematics. The 
latter has two components; the bias and the uncertainty 
(which depends on the curvature of the likelihood func- 
tion, i.e. to what extent does the information matrix con- 
strain the cosmological model in question). As in O'Dea, 
Challinor & Johnson [6] we define 
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I A? 



(26) 



where the superscript refers to values evaluated at the 
peak of the likelihood function, i.e. the values we as- 
sume for the underlying model, and AA^ and Aa\ i are 
the bias (defined in Eq. 24) and the change in the statis- 
tical error for a given experiment and for the parameters 
Xi induced by the beam systematics, respectively. As 
shown in O'Dea, Challinor & Johnson [6] these two pa- 
rameters depend solely on the primordial, lensing and 
systematics power spectra. We require both 5 and (3 not 
to exceed 10% of the uncertainty without systematics. As 
illustrated in Fig. 4 for the case of tensor-to-scalar ratio 
and neutrino total mass, the bias exceeds the uncertainty 
at some value of the beam ellipticity. This is a general 
result; for given beam systematic and a cosmological pa- 
rameter the bias becomes the dominant component of the 
error in parameter estimation for sufficiently large beam 
imperfection (ellipticity, gain, etc). This sets the limit 
on our five systematics parameters as demonstrated in 
Tables IV, V, VI, VII, VIII and IX for PLANCK, PO- 
LARBEAR, SPIDER, QUIET+CLOVER, CMBPOL-A 
and CMBPOL-B (we considered two cases which we re- 
fer to as CMBPOL-A and CMBPOL-B, the former is a 
high sensitivity experiment with 1000 Planck-equivalent 
detectors, the later is motivated by Kaplinghat, Knox 
& Song [12]) whose specifications arc given in Table III. 
For the systematics power spectra we used the expres- 
sions in Table II assuming only temperature leakage (i.e. 



polarization-free underlying sky) except for the effect of 
differential rotation where we consider mixing between 
the E- and B-mode. C TE cannot leak to the B-mode 
power spectrum since it assumes negative values for cer- 
tain multipole numbers and C s -contribution is higher 
order correction to the B-mode systematics and will add 
only few percent at most to the induced systematics. Due 
to the scaling of the systematics with the beam width, 
this potentially negligible C E contribution, which con- 
taminate the B-mode polarization at second order, will 
result in only a few percent change to the tolerance levels 
for the beam parameters we consider. More specifically, 
Tables IV, V, VI, VII, VIII and IX contain the max- 
imum differential gain, beamwidth, pointing, ellipticity 

and rotation, in units of \J~^ of a percent, \J~^ of an 
arcsecond, a percent and a degree, respectively. As ex- 
plained in Eq.(12) the factor 2ir in the denominators is 
the result of our assumption that the scanning strategy 
is spatially uniform. All results are robust against chang- 
ing the step size used. Under the assumption of uniform 
scanning strategy and assuming 'worst case scenario' that 
the quadrupole moment of the scanning strategy is maxi- 
mal {\J~^ = 1) in case of differential gain and beamwidth 
and that the monopole and octupole moments are maxi- 
mal {\J~1^ = 1) these allowed values correspond directly 
to g and \i in (percent units) and p (in arcsec units), re- 
spectively. It is also apparent from the tables that better 
sensitivity and higher resolution experiments generally 
require better control of beam parameters. This is ex- 
pected since our criterion is the fixed 10% threshold in [3 
and 5, i.e. the systematics are not allowed to exceed the 
10% (or any other reasonably chosen threshold) level in 
the parameter uncertainty (<j\ i ) units. Higher sensitivity 
experiments will have smaller <j\ i in general and there- 
fore the allowed g, fi, e, p and e will be smaller. This does 
not imply necessarily that controlling beam systematics 
of higher sensitivity experiments will be more challeng- 
ing since the uncertainty of beam parameters is a direct 
result of the S/N level with which the beam is calibrated 
against a point-source. Reducing the detector noise (akin 
to higher-sensitivity experiment) allows smaller uncer- 
tainty in beam parameters. Also, as mentioned above 
the 10% threshold adapted here is arbitrarily chosen and 
as long as we keep the systematic bias on the cosmolog- 
ical parameters Xi smaller than <j\ i , e.g. even in case 8 
is as high as 0.2, the beam systematics will not signifi- 
cantly degrade the science. Therefore, even very sensitive 
high-resolution experiments are expected to yield good 
systematics control. 
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MCMC Results 

This work is the first to employ Monte-Carlo simula- 
tions to assess the effect of beam systematics on param- 
eter estimation. We first considered beam parameters 
which bias the tensor-to-scalar ratio by 10%, 50% and 
100% of the error (i.e. 6 = 0.1, 0.5 and 1.0, see Eq. 26) 
found with the Fisher matrix formalism. The limiting 
value of S = 0.1 was chosen for the Fisher matrix for- 
malism so that it was small enough such that we should 
not be able to see a bias, which is what we want for a 
limit. For the MCMC simulations, we want to be able 
to distinctly see a bias and with 5 = 1 we expected to 
observe it. The other two values S = 0.1 and 0.5 are also 
reported mainly for the purpose of comparison between 
the naive expectations (based on Fisher matrix analysis) 
with the MCMC results. For these simulations we focus 
only on the most sensitive cosmological parameters: r, w 
and M v . We ran the MCMC analysis on the experiments 
POLARBEAR, CMBPOL-B, and QUIET+CLOVER. 

Our results for POLARBEAR are reported in Table X, 
for CMBPOL-B in Table XI, and for QUIET+CLOVER 
in Table XII. The most sensitive parameter to the beam 
systematics considered here turns out to be r, the tensor- 
to-scalar ratio as it is not constrained by the addition 
of the lensing power spectra. We found the bias on 
r can be as high as ~ 100%. w and M„, which arc 
mostly constrained by the larger (compared to the pri- 
mordial signal from inflation) lensing signal, are changed 
by no more that ~ 10%, itself a non-negligible bias. 
Even in the absence of systematics POLARBEAR and 
QUIET+CLOVER exhibit a small bias (approximately 
ifr) in w (Tables X and XII, respectively) towards values 
smaller than -1 but this situation significantly improves 
with CMBPOL-B. The reason is that, as is evident from 
our simulation, the 1-D distribution for w, while peaked 
at — 1, is skewed towards more negative values. As the 
experiment sensitivity improves, such as in CMBPOL- 
B, this small bias becomes insignificant. Most impor- 
tantly, we also found that the levels of bias (in r) caused 
by differential beamwidth and ellipticity exceed the bias 
found with the naive inclusion of the power spectrum 
bias in the Fisher-matrix formalism (Table XIII). This 
illustrates that the simplistic approach to bias within the 
Fisher-matrix formalism underestimates the induced bias 
on the cosmological parameters. However, the Fisher ma- 
trix can be used, as done here, to determine the starting 
values for MCMC simulations. These two systematics, 
the differential beamwidth and ellipticity, are second or- 
der gradients of the underlying temperature anisotropy 
as opposed to the first order gradient in case of, e.g. first 
order pointing error effect. This implies that for given /i 
and e this effect steeply increases towards smaller scales. 
The Fisher-matrix bias calculation is based, however, on 
the assumption that the bias is relatively small, an as- 



sumption which certainly breaks down when high reso- 
lution experiments are considered (i.e. with SPIDER's 
comparatively low angular resolution, for example, we 
expect the tension between the Fisher-matrix-based and 
MCMC estimations of the bias to be smaller). 

CONCLUSIONS 

The purpose of this work was to illustrate the ef- 
fect of beam systematics on parameter extraction from 
CMB observations. Beam systematics arc expected to be 
significant especially for detecting the B-mode polariza- 
tion. Ongoing and future experiments must meet very 
challenging requirements at the experiment design and 
data analysis phases to assure polarimctric fidelity Ulti- 
mately, a major target of these experiments is the most 
accurate estimation of cosmological parameters, and for 
this end it is mandatory to assess, among other issues, 
the propagation of beam systematics to parameter es- 
timation. The tolerance levels chosen in this work are 
somewhat arbitrary and may be changed at will, accord- 
ing to the goals of individual experiments, and the nu- 
merical values we quote in the tables should be viewed 
in this perspective. 

The only similar work so far to set tolerance levels on 
beam mismatch in the context of parameter estimation 
is O'Dea, Challinor & Johnson [6] which influenced our 
present work. However, we expand on this work in several 
ways. While O'Dea, Challinor & Johnson [6] considered 
only the effect of systematics on the tensor-to-scalar ratio 
r, we consider a family of parameters associated with the 
B-mode sector: r, M v , a, w and f^. We set all other cos- 
mological parameters to be consistent with the WMAP 
values. In order to exhaust the potential of the CMB to 
constrain these parameters we carried out lensing extrac- 
tion. In addition, we repeated the analysis for POLAR- 
BEAR, CMBPOL-B and QUIET+CLOVER with Monte 
Carlo simulations and found that the Fisher-Matrix ap- 
proximation is, in general, inadequate for appraising the 
biases. We also found that high resolution experiments, 
such as POLARBEAR are very sensitive to bias from sec- 
ond order gradient effects (i.e. differential ellipticity and 
differential beamwidth) which is underestimated by the 
Fisher-matrix-based calculation, but fully accounted for 
with MCMC simulations. Also, unlike O'Dea, Challinor 
& Johnson [6] our results are presented independently 
of the scanning strategy details. The only assumption 
we made was that the scanning strategy is spatially uni- 
form, a condition which can be achieved with or without 
a HWP which samples the polarization angles in a way 
which is uniform; both spatially, and in terms of polar- 
ization angle. In case that this approximation fails the 
more general formalism (Shimon et al. [7]) should be 
used with the added complexity introduced to lensing re- 
construction by the scanning-induced non-gaussianity of 
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the systematic B-mode. 

We find that parameter bias is the dominant factor 
and its level actually sets the upper bounds on the beam 
parameters appearing in Tables IV through IX. Our re- 
sults show that the most severe constraints are set on the 
most sensitive experiments for a given tolerance on S and 
(5 since these quantities are experiment-dependent (Eq. 
26) and since, in general, an experiment with higher res- 
olution and better sensitivity will result in smaller errors 
a\ i • We expect that the constraints on the systematics 
should be more demanding so as to realize the potential 
of experiments. As mentioned above, the most stringent 
constraints are obtained from the requirement on the 
bias rather than from increased parameter uncertainty. 
Again, for the same reason, as shown for specific exam- 
ples in Fig. 4; the bias always exceeds the uncertainty 
for large enough systematics and this always takes place 
before the 10% thresholds in Eqs.(26) are attained. The 
reason is that for large enough systematics the induced 
spurious polarization becomes comparable to, or exceeds, 
the underlying polarization signals, therefore biasing the 
deduced value. It is easy to visualize configurations in 
which the bias increases without bound while the 'curva- 
ture' of the likelihood function (i.e. the statistical error) 
with respect to specific cosmological parameters does not 
change. It is also clear from the tables that, in general, 
the tensor-to-scalar ratio is the most sensitive parame- 
ter, and the second most sensitive is a, the running of 
the scalar index (although there are some exceptions). 
If the tensor-to-scalar ratio is larger than the case we 
studied (r = 0.01), this conclusion may change since r is 
mainly affected by the overwhelming B-mode systematics 
on degree scales, a is predicted to vanish by the simplest 
models of inflation and was added to parameter space to 
better fit the WMAP and other cosmological data. As 
is well-known, information from Ly-a systems and other 
LSS probes can, in principle, better constrain a if their 
associated systematics can be controlled to a sufficiently 
accurate level. For these small scales the CMB is not the 
ideal tool to extract information and the error that beam 
systematics induce on a are not significant. 

The upper limits we obtained in this work on the al- 
lowed range of beam mismatch parameters for given ex- 
periments and given arbitrarily-set tolerance levels on the 
parameter bias and uncertainty, constitute very conser- 
vative limits. It can certainly be the case that some of 
the systematics studied here may be fully or partially re- 
moved. This includes, in particular, the first order point- 
ing error which couples to the dipole moment of non-ideal 
scanning strategies (see Shimon et al. [7]). By remov- 
ing this dipole during data analysis the effect due to the 
systematic first order pointing error (dipole) drops dra- 
matically. We made no attempt to remove or minimize 
these effects in this work. Our results highlight the need 
for scan mitigation techniques because the coupling of 
several beam systematics to non-ideal scanning strate- 



gies results in systematic errors. This potential solu- 
tion reduces systematics, which ultimately propagate to 
parameter estimation, and affect mainly the parameters 
considered in this work. A brute-force strategy to idealize 
the data could be to remove data points that contribute 
to higher-than-the-monopole moments in the scanning 
strategy. This would effectively make the scanning strat- 
egy 'ideal' and alleviate the effect of the a priori most 
pernicious beam systematics. This procedure 'costs' only 
a minor increase in the instrumental noise (due to throw- 
ing out a fraction of the data) but will greatly reduce the 
most pernicious reducible beam systematic, i.e. the first 
order pointing error ('dipole' effect). The lesson is clear: 
the rich treasures of cosmological parameters deducible 
from B-mode data require a combination of high polari- 
metric fidelity and judicious data mining. Both are emi- 
nently feasible upcoming CMB polarization experiments. 
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depends on 
beam substructure 


effect 


parameter 


definition 


No 


gain 


9 


3i -32 


Yes 


monopole 


A* 




Yes 


dipole 


P 


Pi - P2 


Yes 


quadrupole 


e 


&x ~ iy y 


No 


rotation 


e 


§(ei +£2) 



TABLE I: Definitions of the parameters associated with the systematic effects. Subscripts 1 and 2 refer to the first and second 
polarized beams of the dual beam polarization assumed in this work. 



effect 


parameter 


ACf E 


ACf 


ACf 


gain 


g 





g 2 h*cf 


g 2 h*cT 


monopole 


A* 







4,x 2 (/a) 4 Cf */i 


pointing 


P 


-c e Jl{lp)Cf '* f 3 


J?(i P )cT*f 2 


Jl{lp)Cf*h 


quadrupole 


e 


-I (z)h(z) C4l Cf 




Hitler 


rotation 


e 





4e 2 Cf 


4e 2 Cf 



TABLE II: The scaling laws for the systematic effects to the power spectra Cf , Cf E , Cf and Cf assuming the underlying 
sky is not polarized (except for the rotation signal where we assume the E, and B-mode signals are present) and a general, not 
necessarily ideal or uniform, scanning strategy. The next order contribution ( 10% of the 'pure' temperature leakage shown 
in the table) is contributed by Cj B . It can be easily calculated based on the general expressions in Shimon et al. [7] where 
the definitions of z, p, e, etc., are also found. For the pointing error we found that the 'irreducible' contribution to B-mode 
contamination, arising from a second order effect, is extremely small and therefore only the first order terms (which vanish in 
ideal scanning strategy) are shown. The functions /1 and fc are experiment-specific and encapsulate the information about the 
scanning strategy which couples to the beam mismatch parameters to generate spurious polarization. In general, the functions 
fi and /2 are spatially-anisotropic but for simplicity, and to obtain a first-order approximation, we consider them constants 
(see sec. 2.2) in general. In the case of ideal scanning strategy they identically vanish. The exact expressions are given in 
Shimon et al. [7] . 
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Experiment 


/sky 


V 


6 b 


At 


A E 


PLANCK 


0.65 


30 


33' 


2.0 


2.8 






44 


24' 


2.7 


3.9 






70 


14' 


4.7 


6.7 






100 


9.5' 


2.5 


4.0 






143 


7.1' 


2.2 


4.2 






217 


5.0' 


4.8 


9.8 






353 


5.0' 


14.7 


29.8 






545 


5.0' 


147 


oo 






857 


5.0' 


6700 


oo 


POLARBEAR 


0.03 


90 


6.7' 


1.13 


1.6 






150 


4.0' 


1.70 


2.4 






220 


2.7' 


8.0 


11.3 


SPIDER 


0.6 


96 


58' 


0.46 


0.65 






145 


40' 


0.50 


0.71 






225 


26' 


2.22 


3.14 






275 


21' 


5.71 


8.08 


QUIET+CLOVER 


0.015 


150 


10' 


0.34 


0.48 


CMBPOL-A 


0.65 


150 


5' 


0.22 


0.32 


CMBPOL-B 


0.65 


150 


3' 


1.0 


1.4 



TABLE III: Instrumental characteristics of the CMB polarimeters considered in this work: / sky is the observed fraction of the 
sky, v is the center frequency of the channels in GHz, 9b is the FWHM (Full- Width at Half-Maximum) in arc-minutes, At is 
the temperature sensitivity per pixel in /jK and A E = A s is the polarization sensitivity. For all experiments, we assumed one 
year of observations (PLANCK [32], POLARBEAR [33], SPIDER [18], QUIET+CLOVER [6]). CMBPOL A & B represent 
toy experiments for illustration with CMBPOL-A having 1000 PLANCK detectors and PLANCK resolution and CMBPOL-B 
has higher resolution but 1/j.K noise level (Kaplinghat, Knox & Song [12]). 



Parameter 


Nominal value 


9 / fl 
1% V 2vr 


M / fl 
1% V 2vr 


( p ) i / h 

V 1" > V 2tt 


e 

T% 


e [deg] 


r 


0.01 


0.02 


0.42 


1.5 


0.8 


0.72 


w 


-1 


0.33 


0.38 


2.5 


2.4 


2.86 


0, k 





0.37 


0.44 


3.0 


2.6 


3.72 


a 





0.67 


0.33 


2.2 


2.1 


2.23 


M„[eV] 


0.56 


0.32 


0.38 


2.4 


2.4 


2.58 



TABLE IV: Systematics tolerance for PLANCK: shown are the nominal cosmological parameters we used along with the 
tolerance levels (as defined by the criterion that both 8 and (3, Eq. 26, should not exceed the 10% threshold) set on combinations 
of the quadrupole of the scanning strategy (/i, under the assumption of uniform scanning strategy) and the dimensionless 
differential gain g and differential beamwidth /i. Also is shown the constraint on pointing error weighted by the dipole of the 
scanning strategy (/2) in arcsec units. The tolerance level on ellipticity e is dimensionless (we assumed the worst-case-scenario 
that if) — 45°) and the allowed rotation e is given in angular degrees. Except for the differential beamwidth effect, the most 
severe constraints are obtained from the requirement that r is not biased, g and fi, the parameters representing the differential 
gain and differential beamwidth, are defined in Table I. 



Parameter 


Nominal value 


_2_,/lL 

1% V 2^ 


-if-, /IT 

1% V 2^ 




e 
1% 


e [deg] 


r 


0.01 


0.01 


0.74 


0.5 


1.4 


0.25 


w 


-1 


0.16 


0.38 


1.7 


1.8 


1.26 


fife 





0.18 


0.39 


1.8 


1.8 


2.01 


Q 





0.17 


0.30 


1.2 


1.3 


0.77 


M„[eV] 


0.56 


0.15 


0.42 


1.9 


1.8 


1.06 



TABLE V: Systematics tolerance for POLARBEAR: As in Table IV. 
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Parameter 


Nominal value 


9 , / /l 
1% V 2ir 


M , ffl 
1% V 2ir 


V 1" I V 2tt 


e 


e [deg] 


r 


0.01 


0.03 


0.10 


2.2 


0.19 


0.97 


w 


-1 


0.13 


0.31 


9.2 


0.47 


2.86 


fl k 





0.10 


0.57 


9.9 


1.75 


3.43 


a 





0.19 


0.12 


5.9 


0.55 


6.88 


M„[eV] 


0.56 


0.10 


0.26 


10.9 


0.38 


3.72 



TABLE VI: Systematics tolerance for SPIDER: As in Table IV. 



Parameter 


Nominal value 


_2_,/lL 

1% V 2ir 


1% V 2ir 




e 

1% 


e [deg] 


r 


0.01 


0.009 


0.20 


0.4 


0.4 


0.2 


w 


-1 


0.114 


0.17 


3.2 


0.6 


0.9 


Q k 





0.122 


0.18 


3.5 


0.7 


1.0 


a 





0.148 


0.13 


1.4 


0.4 


0.6 


M„[cV] 


0.56 


0.109 


0.18 


3.1 


0.7 


0.8 



TABLE VII: Systematics tolerance for QUIET+CLOVER: As in Table IV. 



Parameter 


Nominal value 


S-.fE 

1% V 2tt 


Jf_, in. 

1% V 2ir 


(£)>/& 


e 

i% 


e [deg] 


r 


0.01 


0.0016 


0.05 


0.04 


0.10 


0.023 


w 


-1 


0.0259 


0.19 


0.4 


0.28 


0.773 


Q, k 





0.0270 


0.21 


0.4 


0.28 


0.372 


a 





0.0266 


0.08 


0.3 


0.21 


0.123 


M„[eV] 


0.56 


0.0251 


0.18 


0.4 


0.28 


0.401 



TABLE VIII: Systematics tolerance for CMBPOL-A: As in Table IV. 



Parameter 


Nominal value 


_2_,/lL 

1% V 2ir 


1% V 2?r 




e 

i% 


e [deg] 


r 


0.01 


0.0031 


0.57 


0.2 


l.l 


0.066 


w 


-1 


0.0728 


0.40 


0.9 


1.7 


0.716 


fl k 





0.0762 


0.39 


0.8 


1.8 


0.888 


a 





0.0600 


0.30 


0.4 


1.3 


0.315 


M„[cV] 


0.56 


0.0700 


0.45 


1.4 


1.7 


0.544 



TABLE IX: Systematics tolerance for CMBPOL-B: As in Table IV. 



Parameter 


no sys. 


eFish 


diff. gain 


diff. pointing 


diff. beamwidth 


diff. ellipticity 


diff. rotation 


r 


0.0103 ± 0.0036 


0.1 


0.0108 ±0.0037 


0.0108 ± 0.0038 


0.0109 ± 0.0036 


0.0107 ±0.0038 


0.0109 ± 0.0037 






0.5 


0.0123 ±0.0038 


0.0124 ± 0.0040 


0.0133 ± 0.0042 


0.0126 ±0.0041 


0.0122 ± 0.0049 






1.0 


0.0152 ±0.0042 


0.0152 ± 0.0047 


0.0205 ± 0.0058 


0.0192 ± 0.0056 


0.0147 ± 0.0043 


w 


-1.170 ±0.328 


0.1 


-1.170 ±0.328 


-1.168 ±0.325 


-1.217 ± 0.328 


-1.154 ±0.334 


-1.180 ±0.335 






0.5 


-1.163 ±0.307 


-1.186 ±0.331 


-1.391 ±0.269 


-1.174 ±0.329 


-1.172 ±0.327 






1.0 


-1.171 ±0.321 


-1.194 ±0.334 


-1.575 ±0.217 


-1.153 ±0.333 


-1.145 ±0.320 


M„[eV] 


0.537 ±0.071 


0.1 


0.539 ± 0.068 


0.538 ± 0.070 


0.558 ±0.068 


0.531 ±0.069 


0.540 ± 0.070 






0.5 


0.540 ± 0.067 


0.538 ± 0.073 


0.661 ±0.065 


0.538 ± 0.070 


0.537 ± 0.068 






1.0 


0.544 ± 0.066 


0.542 ± 0.074 


0.924 ± 0.074 


0.533 ± 0.077 


0.533 ± 0.070 



TABLE X: The effect of differential gain, pointing, beamwidth, ellipticity and rotation on parameter estimation for POLAR- 
BEAR obtained with MCMC simulations. The systematic beam parameters e, g, etc. were chosen so that S F " h (Eq.26) 
assumes the specified values (third column from left), i.e. the bias -to-uncertainty ratio in the tensor-to-scalar ratio r (assuming 
r =0.01), as obtained by the Fisher-matrix-based calculation. The values shown are the cosmological parameters recovered 
from the full likelihood function and their la errors. The biases we obtain for differential beamwidth and ellipticity are orders 
of magnitude larger and are not shown here. 



18 



Parameter 


no sys. 


cFish 


diff. gain 


diff. pointing 


diff. bcamwidth 


diff. ellipticity 


diff. rotation 


r 


0.00961 ± 0.00039 


0.1 


0.00964 ± 0.0040 


0.00963 ± 0.0040 


0.00964 ± 0.0039 


0.00962 ± 0.0039 


0.00963 ± 0.0040 






0.5 


0.00982 ± 0.0040 


0.00979 ± 0.00041 


0.00983 ± 0.00040 


0.00980 ± 0.00041 


0.00980 ± 0.00040 






1.0 


0.01004 ± 0.00041 


0.01000 ± 0.00042 


0.01004 ± 0.00041 


0.01002 ± 0.00040 


0.00990 ± 0.00040 


w 


-1.025 ±0.097 


0.1 


-1.027 ±0.095 


-1.027 ±0.097 


-1.045 ±0.100 


-1.025 ±0.095 


-1.029 ±0.097 






0.5 


-1.027 ±0.095 


-1.028 ±0.094 


-1.152 ±0.121 


-1.021 ±0.091 


-1.032 ±0.098 






1.0 


-1.026 ±0.094 


-1.027 ±0.097 


-1.304 ±0.149 


-1.026 ±0.096 


-1.026 ±0.093 


M„[eV] 


0.534 ±0.014 


0.1 


0.534 ±0.014 


0.534 ±0.014 


0.538 ±0.014 


0.534 ±0.014 


0.535 ±0.014 






0.5 


0.535 ±0.014 


0.535 ±0.014 


0.556 ±0.013 


0.534 ±0.014 


0.535 ±0.014 






1.0 


0.535 ±0.014 


0.536 ±0.014 


0.577 ±0.013 


0.534 ±0.014 


0.535 ±0.014 



TABLE XI: The effect of beam systematics on parameter estimation from CMBPOL-B obtained with using MCMC simulations. 



Parameter 


no sys. 


?Fish 


diff. gain 


diff. pointing 


diff. beamwidth 


diff. ellipticity 


diff. rotation 


r 


0.01035 ± 0.00333 


0.1 


0.01076 ± 0.03221 


0.01076 ± 0.00330 


0.01084 ± 0.00323 


0.01053 ± 0.00342 


0.01046 ± 0.00301 






0.5 


0.01199 ±0.00335 


0.01241 ± 0.00365 


0.01287 ± 0.00384 


0.01289 ± 0.00361 


0.01226 ± 0.00361 






1.0 


0.01445 ± 0.00358 


0.01506 ± 0.00409 


0.02326 ± 0.00639 


0.02135 ± 0.00587 


0.01448 ± 0.00386 


w 


-1.143 ±0.362 


0.1 


-1.163 ±0.349 


-1.141 ±0.352 


-1.167 ±0.361 


-1.131 ±0.370 


-1.143 ±0.368 






0.5 


-1.123 ±0.361 


-1.109 ±0.363 


-1.228 ±0.326 


-1.131 ±0.362 


-1.121 ±0.357 






1.0 


-1.137 ±0.356 


-1.153 ±0.362 


-1.347 ±0.322 


-1.160 ± 0.364 


-1.137 ±0.362 


M„[eV] 


0.535 ±0.109 


0.1 


0.536 ±0.105 


0.531 ±0.114 


0.548 ±0.111 


0.522 ±0.112 


0.532 ±0.110 






0.5 


0.526 ± 0.120 


0.526 ±0.114 


0.611 ±0.110 


0.533 ± 0.124 


0.530 ± 0.112 






1.0 


0.528 ±0.117 


0.552 ±0.114 


0.774 ±0.120 


0.530 ±0.123 


0.532 ±0.112 



TABLE XII: The effect of beam systematics on parameter estimation from QUIET±CLOVER obtained with using MCMC 
simulations. 
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Experiment 


Beam parameter 


5 r (5f !sh = 0.1) 


S r (Sp sh = 0.5) 


S r (8^ 3h = 1.0) 


POLARBEAR 


y i 


1 fi 


0.14 


0.56 


1.36 




D\ 
r 


[77 


0.14 


0.58 


1.36 






Pfl 


0.17 


0.83 


2.83 






e 


n 1 1 




9 A7 






£ 


0.17 


0.53 


1.23 


CMBPOL-B 


y 


\l 2ir 


0.08 


0.54 


1.10 




P\ 

r 


[77 
\J 2tt 


0.05 


0.46 


1.15 




pi 


177 


n no 
O.OO 


n re 
U.06 


i in 
1.10 






e 


u.uo 


fl 4Q 


1 (TS 






£ 


0.05 


0.49 


0.74 


QUIET+CLOVER 


y 


^ 2tt 


0.12 


0.49 


1.23 






[17 
\l 2tt 


0.12 


0.62 


1.41 






[77 

V 2tt 


0.15 


0.76 


3.88 






e 


0.06 


0.76 


3.30 






£ 


0.04 


0.57 


1.24 



TABLE XIII: The bias in the tensor-to-scalar ratio r (S r ) obtained with MCMC for POLARBEAR, CMBPOL-B and 
QUIET+CLOVER. These S r values were obtained by assuming each of the five systematics we considered have the values 
which induce a bias S r 7 ' tsh = 0.1, 0.5 and 1, respectively, in the Fisher matrix analysis. This table is a compilation of the cor- 
responding values for r reported in Tables X, XI and XII. Note the discrepancy between the Fisher-matrix-based and MCMC 
forecast for the bias for the differential beamwidth and ellipticity systematics. Both scale as the second-order gradient of the 
temperature Cf as oc l 4 Cf and as a result of this steep rise of the systematics with scale the systematics soon overwhelm the 
primordial B-mode signal and significantly bias the deduced tensor-to-scalar ratio. The Fisher matrix estimate of the bias is 
only a leading order approximation in case the bias is small; an assumption which evidently does not apply to systematics 
which scale as the second order gradient of the temperature anisotropy. 



